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Abstract 

Periodic-disturbance accommodating control is 
investigated for asymptotic momentum manage- 
ment of control moment gyros used as primary ac- 
tuating devices for the Space Station. TLhe pro- 
posed controller utilizes the concepts of quaternion 
feedback control and periodic-disturbance accom- 
modation to achieve oscillations about the constant 
torque equilibrium attitude, while minimizing the 
control effort required. Three-axis coupled equa- 
tions of motion, written in terms of quaternions, 
are derived for roll/yaw controller design and sta- 
bility analysis. The quaternion feedback controller 
designed using the linear-quadratic regulator syn- 
thesis technique is shown to be robust for a wide 
range of pitch angles. It is also shown that the pro 
posed controller tunes the open-loop unstable ve- 
hicle to a stable oscillatory motion which minimizes 
the control effort needed for steady-state operations. 


Introduction 

The Space Station will employ CMGs (control 
moment gyros) as primary actuating devices dur- 
ing normal flight mode operation. Gravity-gradient 
torques will be used for CMG momentum manage- 
ment (unloading). The effect of a constant aerody- 
namic torque on a gravitationally stabilized space- 
craft was first studied by Garber [1]. Such a torque 
produces a constant attitude angle for which aerody- 
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namic and gravitational torques are balanced. Gar- 
ber [1] has shown that small roll/ yaw librational mo- 
tions are affected by large pitch angles. The aero- 
dynamic disturbance torques acting on the Space 
Station are expected to have constant values plus 
periodic components caused mostly by the effects 
of solar panel rotations and Earth's diurnal bulge. 
As a result, attitude and CMG momentum oscilla- 
tion about the torque equilibrium attitude will oc- 
cur. A recent study [2] demonstrates the usefulness 
of the linear-quadratic-regulator synthesis technique 
and the concept of periodic-disturbance accommo- 
dation in minimizing attitude and/or CMG momen- 
tum oscillations as needed for mission requirements. 

This paper is primarily concerned with attitude 
control and CMG periodic-disturbance rejection for 
large-angle pitch maneuvers of the Space Station. 
New results are presented expanding on the con- 
trol scheme developed in [2]. Pitch-coupled roll/yaw 
equations of motion, first discussed in [1], and writ- 
ten in terms of Euler angles, are derived here in 
terms of quaternions. It is shown that these equa- 
tions are well suited for use in designing a roll/yaw 
controller for large pitch motions of the Space Sta- 
tion. A simple concept of using quaternions for the 
control of spacecraft large-angle maneuvers has been 
developed in [3, 4]. The concept is extended here 
to a more complicated case of controlling both the 
attitude and CMG momentum of the Space Sta- 
tion. Furthermore, this paper presents a new con- 
trol concept of asymptotic momentum management 
of the CMGs, which tunes the open-loop unstable 
vehicle to a stable oscillatory motion during steady- 
state operations, while minimizing the control effort 
needed. 

Figure 1 is a functional block diagram represen- 
tation of a quaternion feedback control system pro- 
posed for the Space Station. The attitude determi- 
nation system utilizes rate gyros and star trackers to 
compute inertial quaternions and the absolute angu- 
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lar velocity of the Space Station. Relative quater- 
nions with respect to a local reference frame are then 
computed for control purposes. The proposed atti- 
tude/momentum controller utilizes relative quater- 
nions, body rates (or relative quaternion rates), and 
CMG momenta to generate proper control torque 
commands to the CMGs. 


Space Station Dynamics: 
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Mathematical Models 

In this section, equations of motion for the Space 
Station in a circular orbit are derived in terms of 
quaternions. For simplicity, the Space Station is as- 
sumed to be a single rigid body. Emphasis is on the 
use of quaternions in the equations of motion and in 
feedback control. Quaternions define the rigid body 
attitude as an Euler-axis rotation. The vector part 
of the quaternions indicates the direction of the Eu- 
ler axis. The scalar part of the quaternions is related 
to the rotation angle about the Euler axis. Detailed 
discussion of the kinematics associated with quater- 
nions and Euler angles can be found in [ 5 , 6 ]. 

The relationships between quaternions and Euler 
angles, for the pitch-yaw-roll body-axis rotation se- 
quence used in this paper, are 
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where C\ = cos(0,/2), 5,- = sin(0,/2) for * = 1,2,3, 
(#1, 0 2 , # 3 ) are the roll, pitch, and yaw Euler angles 
of the body axes with respect to the local vertical 
and local horizontal (LVLH) axes, which rotate with 
the orbital angular velocity; and (q u q 2l q 3 ) are the 
vector parts of the quaternions which indicate the 
direction of the Euler axis while q 4 is the scalar part 
of the quaternions and is related to the angle of ro- 
tation about the Euler axis. Inverse relations may 
also be written as: 
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Attitude Kinematics (with respect to LVLH): 
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CMG Momentum: 
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The nonlinear equations of motion and attitude 
kinematics for the Space Station are as follows: 
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and (u>i, u> 2 , w 3 ) are the body-axis components of 
the absolute angular velocity of the Space Station; 

(* = j) are the moments of inertia; (,' ^ _,') 
are the products of inertia; (h u h 2 , h 3 ) are the 
body-axis components of the CMG momentum; (tii , 
« 2 , «a) are the body-axis components of the control 
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torque; (w x , w 2 , w 3 ) are the body-axis components 
of the external disturbance torque; and n is the or- 
bital rate of 0.0011 rad/sec. 

When body and control axes are aligned with 

the principal axes of the Space Station (7i 7 n> 

h = 722, 7 3 = 7 33 ), Eqs. (3) become 
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which can be found in [5]. In this paper, however, 
we present a new set of equations determined by 
linearizing the above equations for the case of large 
pitch angles with small roll/yaw attitude changes. 
In this case, Eqs. (4), (5), and (6) can be linearized 
with respect to 91 and 93 ^ follows: 

Space Station Dynamics: 
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Inverse relations for Eqs. (10c) are written as 
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Equations (7) and (8) may be used to derive the 
Space Station equations of motion in terms of 
quaternions, written as follows: 
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The quaternion relations of Eqs. (10) may be used 
to transform Eqs. (12) to the following known form 
[1, 6] involving only Euler angles: 
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A final linearization with respect to small pitch mo- 
tions leads to the following well-known equations of 
motion: 


For the case with large pitch angles and small 
roll/yaw attitude changes, the relationships between 
quaternions and Euler angles can be simplified by 


Space Station Dynamics: 
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Attitude Kinematics (with respect to LVLH): 
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These linearized equations are used in [2] for the 
case of small roll, pitch, and yaw attitude changes. 
In this paper, emphasis is on the use of Eqs. (7) 
for the momentum/attitude control of the Space 
Station having small roll/yaw attitude changes but 
large-angle pitch motions. 

Inertia values for the Phase 1 Space Station, as 
well as assembly flight 3, are listed in Table 1. Ta- 
ble 2 includes expected aerodynamic disturbances 
which are modeled as a bias plus periodic terms in 
the body- fixed control axes: 

w (t) — Bias 4 A n sin(nt 4 <j> x ) 

-I- A 2n sin(2nt -f <p 2 ) 

+ A 3n sin(3n* 4- <j> 3 ) 

-I- A 4n sin(4nf 4 <£4) 

The disturbance torque acting in each axis is deter- 
mined from data generated at NASA Johnson Space 
Center by a nonlinear simulation program. The pro- 
gram simulates translational and rotational motions 
of the Space Station in orbit about an oblate Earth. 
It includes rotating solar panels, time-varying sur- 
face areas, and time-varying center-of-pressure lo- 
cations. A Jacchia-Lineberry atmospheric model is 
used to compute density variations. The dominant 
aerodynamic torque frequencies at n and 2n are 
caused by Earth s diurnal bulge and solar panel ro- 
tation effects, respectively. Actual magnitudes and 
phases of these disturbance torques are assumed un- 
known for control design. 


Control Issues 

Before presenting the pitch and roll/yaw con- 
troller designs, it is important to clarify some issues 
related to the effects of large pitch motions and in- 
ertia value uncertainties on the stability of the con- 
trolled Space Station. A characteristic of momen- 
tum/attitude control using gravity-gradient torque 
is that pitch, roll, and yaw responses will settle down 


to, or oscillate about, a constant torque equilibrium 
attitude (TEA). Primary factors involved in deter- 
mining the constant pitch TEA are the magnitude of 
the bias in the disturbance torque and the numerical 
difference between roll and yaw moments of inertia. 
This can be seen by studying the steady-state form 
ofEq. (14b): 
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where 0 2 is the pitch TEA angle and w is the bias of 
the pitch disturbance torque. Pitch gravity-gradient 
torque is largest when the pitch attitude is 45 de- 
grees. This is predicted in Eq. (13b) where, at the 
steady-state, 
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It may be necessary to consider this worst-case pitch 
TEA in control system design. After switching to 
CMG mode from some other modes (e.g., reboost 
mode utilizing reaction jets), the Space Station must 
be able to achieve TEA in each axis without CMG 
momentum or commanded torque saturations. 

The large motions possible in the pitch axis em- 
phasize the importance of using Eqs. (7), (12), or 
(13) in designing the control system. These equa- 
tions show the dependence of roll/yaw dynamics on 
pitch attitude. In fact, a roll/yaw closed-loop sys- 
tem designed for small pitch angles may become 
unstable at large pitch angles; therefore, roll/yaw 
closed-loop stability must be checked at various 
pitch TEA values. 

Other factors affecting closed-loop stability are 
uncertainties and variations in moments and prod- 
ucts of inertia. In particular, the magnitude of the 
pitch gravity-gradient torque depends on the differ- 
ence between the roll and yaw moments of inertia. If 
these inertias are nearly equal (e.g., see Table 1, as- 
sembly flight 3 data), pitch gravity-gradient torque 
is very small and the pitch TEA (if it exists) is large. 
In addition, very small uncertainties in the moments 
of inertia can cause the system to become unstable. 
A discussion of the importance of checking closed- 
loop system robustness with respect to inertia un- 
certainties, by varying inertia values in an appropri- 
ate direction”, is included in the appendix of this 
paper. 

Under normal operating conditions, the Space 
Station will have rotating solar arrays. This causes 
time-varying (sinusoidal) roll and yaw moments of 
inertia, and consequently, a similarly time-varying 
gravity-gradient torque in the pitch axis as shown 
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in Fig. 2. If the bias value of Ii(t) - h(t) is such 
that the pitch gravity-gradient torque never changes 
sign, the system will remain stable. If the gravity- 
gradient torque does switch signs, the system may 
become unstable, depending upon how long the sign 
of the torque is changed. For sufficiently short peri- 
ods of this opposite torque, the system will be sta- 
ble with large, bounded responses about the TEA. 
If gravitv-gradient torque is to ultimately be used 
in momentum/attitude control, the above factors 
will be useful in the design of the inertia configu- 
ration. Consideration should be given to defining a 
boundary near the point or condition corresponding 
to zero pitch-axis gravity-gradient torque, defined 
by jj(t) - I 3 (t) = 0. This is illustrated in Fig. 2. 
The boundary may be thought of as the minimum 
gravity-gradient torque allowed in the pitch axis for 
which closed-loop stability, with respect to inertia 
uncertainties, is maintained. 

In the next sections, the pitch and roll/yaw con- 
troller designs are presented along with time sim- 
ulations of the closed-loop system. In this paper, 
we expand on the previous study [2] by developing 
a periodic-disturbance accommodating controller 
which achieves asymptotic momentum management 
of the CMGs in all three axes. The Phase 1 iner- 
tia configuration listed in Table 1 is used. Effects 
of products of inertia are assumed negligible. The 
corresponding time-varying gravity-gradient torque 
in the pitch axis does not change sign; however, the 
time simulations presented here use the assumption 
that inertia values remain constant while periodic 
terms in the aerodynamic torque include solar panel 
rotation effects. Large pitch TEA responses are pro 
duced by introducing an appropriately large bias in 
the pitch-axis disturbance torque. 


Pitch Control 

In this section, a pitch-axis controller is developed 
for attitude and CMC momentum control. It is 
shown in [2] that disturbance rejection filters can be 
used to reject either attitude or CMG momentum 
oscillations occurring at the frequencies present in 
the disturbance torques. Since asymptotic momen- 
tum management of the CMGs in all three axes is 
of primary interest in this paper, the disturbance 
rejection filters for the pitch axis have the following 
forms: 

< i *2 + ( rc ) 2 <*2 = fa 
fa + (2 n) 2 p 2 = fa 


72 4 (3n) 2 72 = fa 
7/2 4- (4/}) 2 7/2 — fa 

where initial conditions for the filter states can be 
arbitrarily selected (usually zero initial conditions). 
Use of filters at frequencies ??, 2 n, 3n, and An is in- 
dicated by aerodynamic torque data generated by a 
nonlinear simulation program written for the Space 
Station. The pitch-axis control logic is given by a 
single control input involving twelve states: 

U'2 = ^2*2 0 ^) 

where 

K 2 = a 1 x 12 gain matrix 

A [ q 2 Q2 fa f fa a 2 fa 

X2 ~~ fa fa 72 72 7?2 ? ?2 Y • 

The control task is to find proper gains for this 

twelve-state feedback controller. 

In order to use linear control design methodolo- 
gies, Eq. (Tb) must be linearized for small pitch mo- 
tion. This results in Eq. (14b) t which is used as the 
basis for pitch control analysis and design. Vari- 
ous techniques may be used in selecting the twelve 
gains of Eq. (16). These include linear-quadratic- 
regulator (LQR) synthesis [7] and direct assignment 
of closed-loop eigenvalues using a pole-placement 
technique. Several iterations of any method may 
be required to achieve satisfactory closed-loop per- 
formance and robustness. Note that gains resulting 
from Eqs. (14b) are for the use of 0 2 in state feed- 
back. In order to accommodate q 2 for use in feed- 
back, the gains corresponding to states q 2 and q 2 are 
doubled since the approximation used for q 2 is 0 2 /2. 
New gains do not need to be computed for the pitch 
controller in the case of large pitch motion. It is 
mostly roll/yaw destabilization at large pitch angles 
which forces pitch-axis instability. 

The open-loop pitch axis of the Phase 1 inertia 
configuration is unstable, with poles at s = ±1.5n, 
0, 0, and filter poles at s = ±jn, ±j2n, ±j37i, ±j47i. 
One pole at s = 0 comes from the integral feedback 
of /i 2 * After iterative use of an LQR synthesis code, 
available in CTRL-C software, a set of closed-loop 
eigenvalues have been selected and are listed in Ta- 
ble 3. The corresponding gain set is given in Table 4. 
Closed-loop pitch responses of Eq. (6), for a pitch- 
axis maneuver of —30 degrees (caused by a large 
pitch-axis torque bias), can be seen in Figs. 4, 5, 
and 6. Comments on the responses are reserved un- 
til after the presentation of the roll/yaw controller 
design. 
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Roll/Yaw Control 


where 


The roll/yaw controller has a structure similar to 
that of the pitch controller. By examining the open- 
loop transfer function matrix from control inputs to 
roll/yaw attitudes and CMG momentum, it is shown 
in [2] that a periodic disturbance at the orbital rate 
can be rejected in the yaw attitude but not in the 
roll attitude. The analysis is accomplished using 
Eqs. (9) and (14), which assume small motions in 
all axes. In an effort to determine if it is possible to 
have periodic-disturbance rejection in both roll and 
yaw attitudes for this different case involving large 
pitch motions, a similar analysis is considered here 
using the pitch-coupled roll/yaw dynamics described 
by Eqs. (13). 

By combining Eqs. (9) and (13), the transfer func- 
tion matrix from (ui,u 3 ) to (#i,0 3 ) can be written 
as 

[ ^i( s ) 1 [ Gi^s) G i3 (s) 1 T tii(s) 

L ^3(s) J A [ G 31 (s) G 33 (s) J [ u 3 (s) 

where 

G u (s) = - [/ 3 s 2 + (1 + 3si)n 2 (/ 2 - /j)] 
x [s 2 + n 2 ] 

Gi 3 (s) = - [n(Ji - I 2 + I 3 )s 

- 3 n 2 (7 2 - / 3 )s 2 c 2 ][s 2 + n 2 ] 

G 3i(«) = [n(/i - / 2 + I 3 )s 

+ 3n 2 (7 2 — / 1 )s 2 c 2 ][s 2 -f n 2 ] 

G 33 (s) = - [7 x s 2 + (1 + 3c|)n 2 (7 2 - 7 3 )] 
x [s 2 + n 2 ] 

A = / 1 / 3 (s 2 +n 2 ){s 4 + n 2 [l + 3Jb 1 +ife 1 fc3 

+ 3(k 3 -k 1 )s 2 2 ]s 2 

+ [3n 3 (fc 3 - Ai)s 2 c 2 ]s + 4n 4 Jfci(b 3 } 

and s 2 = sin 0 2 , c 2 = cos 0 2 ,ii = ( I 2 - h)fh, and 
^3 = (^2 — h)/ hi ■ Transmission zeros at ±jn ap- 
pear in the transfer function matrix. It would seem 
that periodic disturbances of frequency n cannot be 
rejected in either the roll or yaw attitudes; however, 
for CMG momentum and control torque relations 
defined by 

hi = uj 

h 3 = u 3 


G n = (s){-I 3 s 3 + n 2 [3(/i - I 2 )s\ - 7 3 ]s 
-I- 3 n 3 (7 2 — 7 3 )s 2 c 2 ) 

G 13 = ( ns ){(f 2 — h)s 2 + [3n(7 2 — 7 3 )s 2 c 2 ]s 
+ (1 + 3« 2 )n 2 (7 2 — 7i)} 

G 31 = - (ns){(7 2 - 7 3 )s 2 + (1 + 3c 2 )n 2 (7 2 - 7 3 )} 

G33 = - {hs 4 + n 2 [3(7 2 - 7 3 )c 2 + 7JS 
+ 3n 2 (7 2 — 7 !)s2C 2 ) 

^ = hil 3 {s ){s 4 + n 2 [l + 3 ki + l*i k 3 

+ 3(k 3 - ki)s 2 ]s 2 

+ [3n 3 (fc 3 - fc 1 )s 2 C 2 ]s + 4n 4 k 1 fc 3 }. 

Transmission zeros are not apparent in these expres- 
sions. A numerical analysis reveals, however, that 
there are transmission zeros at ijn for the trans- 
fer functions from u x and u 3 to roll attitude, while 
yaw attitude has no troublesome zeros. These re- 
sults show that, even for the case of pitch-coupled 
roll/yaw dynamics, there is an inability to reject roll 
attitude oscillations occurring at the orbital rate. 
Hence, in this paper, periodic-disturbance rejection 
for CMG momentum in both the roll and yaw axes is 
considered. That is, a periodic control of the Space 
Station for asymptotic momentum management of 
the CMGs in all three axes is of primary interest 
here. 

Periodic-disturbance rejection filters for the 
roll/yaw axes can be represented as: 

<*1 + {n) 2 ai = li\ 

Pi + (2 n) 2 pi = hi 

Ti + (3n) 2 7! = hi 
iji + (4n) 2 /7, = h 1 
«3 + (n) 2 a 3 = h 3 
P 3 + (2 nfp 3 = h 3 
73 + (3n) 2 7 3 = h 3 
Vs + (4 ti) 2 t ) 3 = h 3 . 

The roll/yaw control logic involving two control 
inputs and twenty-four states is expressed as 



where 


and with appropriate alterations of Eqs. (13), the 
transfer function matrix from (u x ,u 3 ) to (0 X ,0 3 ) can 
be written as 


' Oi(s) = I r GiiW G l3 (s) 

. ^ 3 ( 5 ) A [ £ 31 ( 5 ) G 33 (s) 


«i(«) 

U 3 (s) 


K = a 2 x 24 gain matrix 

A [ <71 ui h\ f hi ai 

Pi Pi 7i 7i m 7)1 ] 

X 3 = ^ ^3 

Pi 3 Pz 73 73 rf 3 773 


T (roll states) 
(yaw states) . 
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The pitch-coupled roll/yaw equations described by 
Eqs. (7) are used as the basis for the roll/yaw con- 
troller design. It is evident in the pitch-coupled 
equations that roll/yaw dynamics are dependent 
upon pitch attitude. Consideration must be given 
to the pitch-axis TEA when designing the roll/yaw 
controller gains. Equations (7) (or Eqs. (12) or (13)) 
are especially useful for this purpose when q 2 
and 74 in the equations are assigned their respec- 
tive values corresponding to the expected pitch 
TEA (see Eqs. (10a) and (10b)). As a result, 
Eqs. (7a) and (7c) become linear and any linear con- 
trol design methodologies may be used to design the 
roll/yaw controller gains. 

For spacecraft operating under the influence of 
gravity-gradient torques, it is interesting to exam- 
ine the changes in roll/yaw open-loop eigenvalues 
that occur as the pitch bias changes. 1 his was first 
studied by Garber [1]. A root locus of open-loop 
eigenvalues versus pitch angle, for the Phase 1 in- 
ertia configuration, is shown in Fig. 3. It can be 
seen that the open-loop roll/yaw dynamics are not 
very sensitive to pitch attitude. The Space Station 
is unstable with poles at s = ±1.05n±j0.7n, 0, 0, 
ijn, and filter poles at s = ±jn, ±jn, 4:j2n, ±j2n, 
±j3n, ±j3n, ±j4n, ±j4n (for 0 2 = 0° where q 2 ~0 
and 94=1). The double pole at s = 0 occur because 
of the integral feedback of h\ and /13. 

After iterative use of an LQR synthesis code, 
closed-loop eigenvalues have been selected and are 
listed in Table. 3. A gain set for Eqs. (17), corre- 
sponding to a pitch TEA of 0°, is listed in Table. 4. 
For these gains, the closed-loop roll/yaw axes are 
stable for pitch angles ranging from —21° to 4-23°. 
Since a simulation of the large-angle pitch maneu- 
ver needed to reach a pitch TEA of —30° is desired, 
a different gain set is used for the simulations pre- 
sented in this paper. For these gains, the closed- 
loop roll/yaw axes are stable for pitch angles ranging 
from —48° to 4-3°. Closed-loop roll/yaw responses 
of Eqs. (6), for a pitch-axis maneuver of —30 de- 
grees, can be seen in Figs. 4, 5, and 6. The overall 
closed-loop system has a 10 dB gain margin and a 
phase margin of 60° in each control loop. 


Discussion of Simulation Results 

Closed-loop responses for a simulation of the non- 
linear dynamics described by Eqs. (6) are shown in 
Figs. 4, 5, and, 6. Quantities plotted include quater- 
nions, CMG momenta, and control torques. Initial 
conditions corresponding to 0i(O) = #2(0) = ^3(0) “ 


1 degree and #i(0) = 0 2 (0) — ^3(0) = 0-001 deg/sec 
are assumed. Allowable limits on CMG momentum 
and commanded torque are assumed to be 30,000 
ft-lb-sec and 150 ft-lb, respectively. 

In the roll axis, quaternion q\ oscillates (±0.7°) 
about a roll TEA of -0.003 (0 l = -0.5°), Roll 
CMG momentum h\ is the input to the roll-axis 
disturbance rejection filter and settles down to zero 
after reaching a maximum value near 8000 ft-lb-sec. 
Control torque uy is zero at the steady-state and 
has a peak value near 23 ft-lb. In the pitch axis, 
quaternion q 2 oscillates (±4.3°) about a pitch TEA 
of -0.257 [O'? “ -30°). The large-angle maneuver 
causes the pitch CMG momentum ho (disturbance 
filter input) to become quite large at nearly 18,000 
ft-lb-sec, before settling to zero. Control torque u 2 is 
zero at the steady-state with a maximum value near 
30 ft-lb. In the yaw axis, quaternion oscillates 
(±1.2°) about a yaw TEA of —0.013 ($3 = —1.5°). 
Yaw CMG momentum /13 (disturbance filter input) 
settles down to zero after reaching a maximum value 
close to 000 ft-lb-sec. Control torque u$ is zero at 
the steady-st ate, and reaches a maximum value near 
12 ft-lb. 

The simulations show that the proposed control 
scheme tunes the open- loop unstable Space Sta- 
tion to a stable, oscillatory motion which mini- 
mizes control effort during steady-state operations. 
For the assumed disturbance torque models (with 
unknown magnitudes and phases), the stabilized 
Space Station needs no control torque at steady- 
state conditions. Analysis shows, however, that 
small-amplitude periodic components of frequencies 
5 71 and 6n are present in 1/1 and 113 at the steady- 
state. These small residual components are caused 
by the coupling between the pitch and roll/yaw axes, 
and become particularly noticeable for large pitch 
biases. 

There is an interesting feature of the quaternion 
feedback scheme which is not apparent from the sim- 
ulation responses. If rejection of pitch and yaw at- 
titude oscillations is desired, it seems natural to use 
q ? and 73 as inputs to the respective disturbance re- 
jection filters. Even though q 2 and q^ will become 
constant at the steady -state, all of the Euler angles 
will oscillate. A study of Eqs. (2) (or Eqs. (11)) re- 
veals why. By assigning constant values to q 2 , q 3, 
and q A in Eqs. (2) (or Eqs. (11)), it can be seen that 
0 y< 0o , and 0 3 are all functions of qy, which oscil- 
lates. The same may he said for q \ , q 2 , 73, and if 
Euler angle feedback is used. For oscillations of 0\ 
and constant values for 0 2 and Eqs. (1) show that 
all quaternions are functions of 9\ and will therefore 
oscillate. In either case, however, these oscillations 
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are small. The important point is that the elimi- 
nation of pitch or yaw oscillations, if needed, may 
be accomplished by using Euler angles 0 2 and 0 3 as 
disturbance rejection filter inputs (with appropriate 
gain changes). 

A check of closed-loop robustness with respect to 
inertia uncertainties emphasizes an important issue 
associated with the Phase 1 inertia configuration. 
By selecting specific “directions’' in which to vary 
the three moments of inertia (I 1} / 2 , / 3 ), the closed- 
loop system can be shown to be unstable for as lit- 
tle as — i % uncertainty in I 3 with T8% uncertainty 
m 1 1 . For these inertia variations, the pitch-axis 
gravity-gradient torque disappears (Ij — / 3 = 0) and 
closed-loop pitch dynamics become unstable. The 
limitations shown in this example (and several oth- 
ers involving even smaller inertia uncertainties), are 
not related to the selection of control logic but are 
physical limitations inherent to the inertia configu- 
ration of the Phase 1 Space Station. A description 
of the inertia variation “directions” used above is 
presented in the appendix of this paper. 


Conclusions 

In this paper asymptotic momentum manage- 
ment of control moment gyros of the Space Sta- 
tion has been investigated. It was shown that the 
proposed controller tunes the Space Station, which 
has a gravitationally “unstable” inertia configura- 
tion, to a stable, oscillatory motion which minimizes 
the control effort needed at the steady-state. By 
utilizing the concepts of quaternion feedback con- 
trol and periodic-disturbance rejection filters, the 
proposed controller provides robust control of the 
Space Station for large-angle pitch motions. The 
pitch-coupled roll/yaw equations of motion derived 
in this paper were shown to be particularly useful in 
roll/vaw controller design and stability analysis. 


Appendix: Inertia Variation Directions 

It is a common practice in control design to satisfy 
time and frequency-domain requirements first, then 
check for closed-loop robustness. For spacecraft, in- 
ertial properties may be very sensitive parameters 
m the closed-loop system. It seems reasonable to 
increase and decrease all inertias by the same per- 
centage, thus checking controller effectiveness for a 
proportionally heavier and lighter spacecraft. This 
procedure may not indicate the true inertia sensi- 


tivity. It is important to consider the magnitude 
and direction of the variation for each inertia value. 
Since the inertia matrix may be transformed to three 
principal moments of inertia by aligning the body 
and principal axes, suggestions for variations in the 
moments of inertia for the roll ( ) , pitch (/ 2 ), and 
yaw (/ 3 ) axes are presented here. 

Three important relationships may be derived 
from the definitions for the moments of inertia. 
These relationships are as follows: 

h T h > 7 3 

h T I 3 > 1 2 
I 2 + I 3 > 1 \ . 

Together, these relations define the physically pos- 
sible inertia configurations. A control designer may 
unknowingly use inertia variations which result in 
inertia values that violate these constraints. Sta- 
bility of the closed-loop system will be tested for 
a physically impossible inertia configuration. The 
important point is not the fictitious inertias but 
whether or not the control designer can redirect this 
extraneous stability margin to encompass more of 
the region of physically possible inertia values. 

When gravity-gradient torque is used in the con- 
trol of a rotating spacecraft, additional inertia con- 
straints are introduced. The control scheme pre- 
sented in this paper is a good example. Equa- 
tions (14) show that roll-axis gravity-gradient and 
gyroscopic coupling torques are zero when I 2 = / 3 , 
pitch-axis gravity-gradient torque is zero when I\ — 
I 3 , and the yaw-axis gyroscopic coupling torque is 
zero when I\ — I 2 . 

A useful aid for visualizing the relationship be- 
tween inertia constraints and inertia variations 
is now presented. Figure A.l shows a three- 
dimensional, cubic figure defined in three “inertial” 
directions. The inertia constraint relations may now 
be visualized as planes in this “inertial” space. The 
planes I x + I 2 = / 3 , l x + / 3 = / 2 , and I 2 + I 3 = J x 
are labeled in Fig. A.l, and represent the physical 
boundaries of inertia values. The area inside the 
three intersecting boundaries represents all physi- 
cally possible inertia configurations. A representa- 
tion of the cut-away portion of the “inertia cube” is 
shown in Fig. A.2. Planes defining the physical and 
system torque boundaries are labeled. Figures A.l 
and A. 2 provide a three-dimensional representation 
of the information presented in the k 3 versus k\ in- 
ertia ratio plots in [5-7]. It may be convenient to 
n ormalize the moments of inertia being studied by 

T /f 4- / 3 in order to locate the position of the 
nominal configuration within a “unit inertia cube”. 
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The relative positioning of the nominal inertia con- 
figuration from the constraint boundaries can then 
be easily determined. 

Since the shortest distance from a point to a plane 
is in a direction perpendicular to that plane, it seems 
logical to check inertia variations in directions per- 
pendicular to the inertia boundaries. In this way, 
the minimum variation necessary to reach a physi- 
cal boundary can be found while checking the closed- 
loop stability of the system in question. For rotating 
spacecraft with gravity-gradient control systems, in- 
ertia variations perpendicular to the planes /i = h, 
Iy = / 3) and I 2 = ^3 are needed. It may be seen 
in Fig. A.2 that these planes intersect inside the re- 
gion of physically possible inertia values, and par- 
tition the region into several sections. It should be 
a control designers’ goal to include the area within 
the physical boundaries inside a “control surface” 
which contains all of the inertia values for which the 
closed-loop system is stable. 
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Table 1. Space station inertia 
configurations 


Inertia 

(slug-ft 2 ) 

Assembly 

Flight#3 

Phase 1 

In 

23.22E6 

50.28E6 

1 22 

1.30E6 

10.80E6 

1 33 

23.23E6 

58.57E6 

1 12 

— 0.023E6 

— 0.39E6 

1 13 

0.477E6 

-0.24E6 

I 23 

— 0.011E6 

0.16E6 


Table 2. Phase 1 aerodynamic torque models 
(in units of ft-lbs) 

w x l+sin(nt)+0.5sin(2nt) 

+0 3sin(3nt)+0.5sin(4nt) 
W 2 13* + 1.2sin(nt)+3.5sin(2nt) 

+0.3sin(3nt)+0.5sin(4nt) 
w 3 l+sin(nt)+0.5sin(2nt) 

+0.3sin(3nt)+0.5sin(4nt) 

* nominal pitch bias torque is 4 but 13 is used 
to produce a large pitch TEA 
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Table 3. Phase 1 closed-loop eigenvalues 
(in units of orbital rate - 0.0011 rad/sec) 



Momentum/ Attitude 

Disturbance Filters 

Pitch 

-1.0, -1.5 

— 1.5±jl.5 

— 0.3ijl.0 
— 0.3±j3.0 

— 0.3±j2.0 
— 0.3±j4.0 

Roll/Yaw 

-0.23,-0.71 
— 1.04±j0.70 

-0.53±jl.54 

-1.06±j0.71 

— 0.14±j0.99 

-0.19±j2.01 

— 0.32±j3.02 
— 0.53±j3.97 

-1.13±j0.75 

-0.47±j2.20 

— 0.6S±j3.21 
— 0.25±j4.00 


Table 4. Phase 1 controller gains 


Pitch 

Units 

Roll/Yaw* 

[K 2 f 


I (K] r 

4.2425E+2 

(ft-lb/rad) 

3.8526E+3 

3.7381E+2 

2.5412E+5 

(ft-lb-sec/rad) 

1.2003E+6 

1.0126E+5 

1.4840E— 2 

(ft-lb/ft-lb-sec) 

1.4360E— 2 

1.9364E— 3 

4.0150E— 6 

(ft-lb/ft-lb-sec 2 ) 

— 1.6361E— 6 

2.7852E— 7 

— 1.9064E— 9 

( ft-lb-rad 2 /ft-lb-sec 3 ) 

3.6578E— 10 

— 1.8526E— 10 

2.1970E— 6 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

7.6282E— 7 

1.1857E— 7 

— 4.6097E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

— 3.2712E— 9 

— 5.7517E— 10 

— 5.2383E— 7 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

— 3.3865E— 7 

— 1.6409E— 7 

— 5.3793E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

— 1.0702E— 8 

— 1.1317E— 9 

— 1.9423E— 6 

(ft-lb-rad 2 / ft-lb-sec 2 ) 

— 3.4827E— 6 

— 5.3664E— 7 

— 7.3458E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

— 1.5903E— 8 

— 1.5491E— 9 

-2.6056E-6 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

— 3.1256E— 6 

— 4.7197E— 7 



9.4016E+2 

2.4994E+2 



— 1.2743E+5 

1.1386E+5 



— 2.4992E— 3 

— 3.5209E— 3 



— 7.3398E— 7 

— 1.0348E— 6 



! 4.8557E— 9 

— 5.5935E— 10 



3.7017E— 7 

— 4.2651E— 6 



2.0608E— 9 

— 6.8224E— 10 



1.8854E— 6 

— 2.4769E— 6 



4.0142E— 10 

9.4962E— 10 



1.5548E— 6 

— 2.7820E— 6 



8.3363E— 10 

8.3453E— 10 



1.3125E— 6 

— 2.5757E— 6 


* Designed for a pitch TEA of 0° 


Units 


(ft -lb/ rad) 
(ft-lb-sec/rad) 

(ft-lb/ ft-lb-sec) 
(ft-lb/ft-lb-sec 2 ) 
(ft-lb-rad 2 / ft-lb-sec 3 ) 
(ft-lb-rad 2 / ft-lb-sec 2 ) 
(ft-lb-rad 2 /ft-lb-sec 3 ) 
(ft-lb-rad 2 /ft-lb-sec 2 ) 
(ft-lb-rad 2 /ft-lb -sec 3 ) 
(ft-lb-rad 2 /ft-lb-sec 2 ) 
(ft lb-rad 2 /ft-lb-sec 3 ) 
(ft-lb-rad 2 /ft-lb -sec 2 ) 
(ft-lb/rad) 
(ft-lb-sec/rad) 
(ft-lb/ft-lb-sec) 
(ft-lb/ft-lb-sec 2 ) 
(ft-lb-rad 2 / ft-lb-sec 3 ) 
(ft-lb-rad 2 /ft-lb -sec 2 ) 
(ft-lb-rad 2 /ft-lb-sec 3 ) 
(ft-lb-rad 2 /ft-lb-sec 2 ) 
(ft-lb-rad 2 / ft-lb-sec 3 ) 
(ft-lb-rad 2 / ft-lb-sec 2 ) 
(ft-lb-rad 2 /ft-lb-sec 3 ) 
(ft-lb-rad 2 / ft-lb-sec 2 ) 
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r 1 



Momentum L j 


B/I indicates body frame with respect to inertial frame 
B/L indicates body frame with respect to LVLH frame 
I/L indicates inertial frame with respect to LVLH frame 

Fig. 1. A quaternion feedback control system 
for the Space Station. 


Crovity-Cradienl Torque Derivative 



Fig. 2. Time- varying pitch-axis gravity-gradient 
torque derivative. 
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Fig. 3. Plot of roll/yaw open-loop poles versus pitch angle 
in the s plane (s in units of orbital rate). 
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Fig. 4. Closed-loop quaternion responses. 
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Fig. 5. CMG momentum histories. 


TIME (ORBITS) 


i U. yy 


Fig. 6. Commanded torque histories. 
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Fig. A.l. Physical boundaries of inertia values. 


Fig. A. 2. Regions within the physical boundaries. 
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